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Abstract 

Exclusive double difFractive events (EDDE) are considered in the framework of the 
Regge-eikonal approach and perturbative calculations for "hard" subprocesses. Total 
and differential cross-sections for processes p + p—^p + X+ p are calculated. 
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1 Introduction 



LHC collaborations aimed at working in low and high regimes related to typical un- 
dulatory (diffractivc) and corpuscular (point-like) behaviours of the corresponding cross- 
sections may offer a very exciting possibility to observe an interplay of both regimes [1]. 
In theory the "hard part" can be (hopefully) treated with perturbative methods whilst 
the "soft" one is definitely nonperturbative. 

Below we give several examples of such an interplay: exclusive particle production by 
diffractively scattered protons, i.e. the processes p + p p + X +p. where -f means also a 
rapidity gap and X represents a particle or a system of particles consisting of or strongly 
coupled to the two-gluon state. 

This process is related to the dominant amplitude of exclusive two-gluon production. 
Driving mechanism of the diffractive processes is the Pomeron. Data on the total cross- 
sections demands unambiguosly for the Pomeron with larger-than-one intercept, thereof 
the need in " unitarisation" . 

As will be seen below, EDDE gives us unique experimental possibilities for particle 
searches and investigations of diffraction itself. This is due to several advantages of the 
process: a) clear signature of the process; b) possibility to use "missing mass method", 
that improve the mass resolution; c) background is strongly suppressed; d) spin-parity 
analysis of the central system can be done; e) interesting measurements concerning the 
interplay between "soft" and "hard" scales are possible. 

2 Calculations 

In Figs. 1,2 we illustrate in detail the process p-\-p^p-\-X-\-p. Off-shell proton-gluon 

amplitudes in Fig. 1 are treated by the method developed in Ref. [2], which is based on 
the extension of Regge-eikonal approach, and succesfuUy used for the description of the 
HERA data [3]. 

The amplitude of the process p + p p-\- X -\-p can be obtained in the following way 
(see Figs. 1,2). The first step is to calculate the "bare" amplitude Tx, which is depicted in 
Fig. 1. The "hard" part is the usual gluon-gluon fusion process calculated by perturbative 
methods in the Standard Model or its extensions. "Soft" amplitudes T12 are obtained in 
the Regge-eikonal approach. The second step is the unitarization procedure, that takes 
into account initial and final state interactions (see Fig. 2). 

We use the following kinematics, which corresponds to the double Regge limit. It is 
convenient to use light-cone components (-I-, — ;-L). The components of momenta of the 
hadrons in Fig. 1 are 
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<i = q) , 

gi = ? + - pi = g + Ai , 
g2 = -g + P2 - P2 = + A2 , 

^1^2 are fractions of protons momenta carried by gluons. For two-dimensional transverse 
vectors we use boldface type. From the above notations we can obtain the relations: 

h,2 = A?,2 - — ' ' ; ' - (2) 



COS 00 — 



-A? 2 , 6,2 ^ 
A1A2 



IA1IIA2 



Mx ^ + tl + t2 - 2A/t]t2 COS 00 

(Pi + of — rri^ + + V2sq^ = Si 
{P2 - — + q^ - \/2sq+ = S2 ■ 

Physical region of diffractive events with two rapidity gaps is defined by the following 
kinematical cuts: 

0.01 GeV^ < \ti,2\ < 1 GeV^ , (3) 
2 

- < 6,2 < Cmax = 0.1 , (4) 



t . ^ 

•,rmn 



si. 



max 



(\/^-v^)'< <(v^+\/^)' (5) 
= iii2S - Ml<^Ml 

We can write the relations in terms of yi,2 and yx (rapidities of hadrons and the system 
X correspondingly). For instance: 

My , 

6,2 c^^e^y-, (6) 

I I ^ 1 / V "^'^ max \ 

\yx\ < 1/0 = In' 



\ Mx ) ' 

yo ~ 2.5 for ^ = 14 TeV , Mx = 0(100 GeV) , 

1 (1 - 6,2)ls 

h,2 



\vi2\ = - In — '■ >9 



In standard terms the amplitude corresponds to the so-called nonfactorized scheme [4] . 
The contribution of the diagram depicted in Fig. 2 is obtained by integrating over all 
internal loop momenta. It was shown in [4], that the leading contribution arises from 
the region of the integration, where momentum q is " Glauber- like" , i.e. of the order 
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(k+m^/\/s, k_m^/y^, km), where k's are of the order 1. The detailed consideration of 
the loop integral hke 

f A /(g,pi,p2, Ai, A2) 
J (27r)4(g' + «0)(g2 + i0)(g2 + i0) 

shows that the main contribution comes from the poles at 

^2 



In this case 



where 



(—^, -^,ql > (8) 



qi2=q2 + Ai^ + 2|q||Ai|cos((/)+^), 
q22=q2 + A2^-2|q||A2|cos(0-^) 
Taking the general form for T-amplitudes that satisfy identities 

= 0, q^T^a = , (9) 
and neglecting terms of the order o(^i), the following expression is found at \ti\ < 1 GeV^: 

K(P, q, qi) = (^^.a - ^) T^,-.,M^ U, qq,) (10) 



^ qqi 

m - _ El 
qqi 

For T^p^gp we use the Regge-eikonal approach [2, 5]. At small U it takes the form of the 
Born approximation, i.e. Regge factor: 

e-^ '^ e^^S (11) 

•So — 55i — J 

where ap(0) = 1.203, «p(0) = 0.094 GeV^^ r^^ = 2.477 GeV^^ are fixed parameters 
for the "hard" Pomeron [5], which have been obtained from the global fit to the data on 
diffractive pp{pp) scattering. Parameters Cgp ~ 3.5, r^^ = 2.54 GeV~^ are defined from 
fitting the HERA data on elastic J/^ production [6], which will be published elsewhere. 
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The upper bound for the constant c^^ ~ 2.3(3.3) can be also estimated from the exclusive 
double diffractive di-jet production at Tevatron (see (41)), if we take CDF cuts and 
the upper limit for the exclusive total di-jet cross-section [7]. The effective value Cgp — 
2.3 corresponds to the case, when the Sudakov suppression factor is absorbed into the 

constant, and Cgp 3.3 IS obtained when taking into account this factor explicitely. 

The full "bare" amplitude looks as follows: 

^ f d'q 8F^-(gi,g2)T/^Jpi, q, gi)Tj^Jp2, -q, ^2) 

ipp^pxp- J ^2^^4 (^2 + ,o)(g2 + ^0)(g2 + ^o) ' ^ ' 



where 



i^n?i,?2) = (^^'^-^)F,,_ 



Factor 8 arises from the colour index contraction. Let = — ~ q^, yx —< yx >— 
and contract all the tensor indices, then the integral (12) takes the form 



T,,^,x, c cle^'''^'^^^, [-^ j ■ 8F,,^x • / , (13) 
b = a'p{0)ln(^]+bo, (14) 



ap(0) 



where Sq — rn? ~ 1 GeV^ is the scale parameter of the model that is used in the global 
fitting of the data on pp{pp) scattering for on-shell amplitudes [5]. It remains fixed in the 
present calculations. It worth nothing that the " rescattering corrections" for the off-shell 
gluon-proton amplitudes, T^, are small (in accordance with a general analysis in Ref. [2]). 
Contrary to this, the "outer" corrections (see Eq.(18)) are significant. 

If we take into account the emission of virtual "soft" gluons, while prohibiting the real 
ones, that could fill rapidity gaps, it results in the Sudakov-like suppression [8]: 



Fs{P) = exp 



3 

'2^ 



(16) 



and in the new value of the integral (15): 

dp 



Jo 



■Fs{P) 



P 



2ap(0) 



^so-m2 + Z72^ 



(17) 



In this case the total cross-section becomes smaller, than without the factor Fg. It plays 
significant role for large Mx- 

Unitarity corrections can be estimated from the elastic pp scattering by the method 
depicted in Fig. 2, where 
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T 



pp~^pXp 



(18) 



V{s , qy) = 4s(27r)252(qy) + 4s 1 rf^be^^Tb 

1 /■ (i^q^ (i^q^ 
iss' J 



V{s , qy) • Tx{pi - qT:P2 + Qt, Ait, Aar) 



less'i (27r)2(27r)2 
• Vis' , q^) , 
Air = Ai-qt-Qt , 
A2T = A2 + + > 

and Spp^pp can be found in Ref. [5]. These "outer" unitarity corrections reduce the 
integrated cross-section by the factor about 14 for the given kinematical cuts and lead to 
the changes in the 0o-dependence. 



3 Results for resonance production 

We have the following expression for the differential cross-section in case of one particle 
production: 



da 7r|T„^™*- 12 



pp—>-pXp I 



dtidt2dCidC2 8s(27r)5V=A 

A ^ k'^ + 2{ti+t2)K+{ti-t2f <0 



(19) 



By partial integrating (19) we obtain t and ^ distributions. The first result of our 
calculations is depicted in the Fig. 3. The antishrinkage of the diffraction peak with 

increasing mass Mx is the direct consequence of the existence of the additional hard scale 
Mx, which makes the interaction radius smaller. The ^ distribution is shown in Fig. 4. 

We can use the following replacement to obtain the cross-section for the EDD process 
p + p—>-p + X+p: 

\Fgg^x\^ ^ AnMxTiX ^ gg) . (20) 
It is possible to simplify calculations after reduction of (19) to 



Tr\T^1'lit. |2 

aa T^\'Lpp^pXp\yx=0 



Ayx , (21) 



dtidt2d(j)o 8s2(27r)5 
where Ayx — 2yo, 0o is the azimuthal angle between outgoing protons. 

4 Standard model Higgs boson production 

For the Standard Model Higgs boson [9] 
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iO 



+ Li", 



.1 + 0" 



1 51 P 



(23) 



where rj — Mfj/Aml, Gp is the Fermi constant, rrit is the top quark mass. NLO K- factor 
1.5 for the gg ^ H process is included to the final answer. 
Numerical results are the following [10] 



Cgp 


Mh (GeV) 




LHC 


TeVatron 


no Sud. suppr. 


Sud. suppr. 


no Sud. suppr. 


Sud. suppr. 


3.5 


100 ^ 500 


110 ^ 57 


4.6 0.14 


12 0.4 


0.5 ^ 0.001 


2.3(3.3) 


100 ^ 500 


20 11 


3.6 ^ 0.11 


2.2 ^ 0.08 


0.4 ^ 0.0009 



We consider four different cases only for Higgs boson as illustration of the Sudakov 
suppressing action. In other examples we take Cgp = 3.3 with Sudakov-like suppression 
according to the CDF data estimations to obtain the lower value of cross-sections at LHC 
and Tevatron. For this case our result is quite close to the one of [8], where the value 
of the total cross-section is about 3 fb. In both cases the most important suppressing 
in the mass region Mh > 100 GeV is due to (perturbative) Sudakov factors, while the 
nonperturbative (absorbtive) factors play relatively minor role. 

Results of other authors were considered in details in [11]. Here we refer to the highest 
cross-section 2 pb for Mh = 400 GeV at LHC energies that was obtained in Ref. [12]. A 
nonfactorized form of the amplitude and a " QCD inspired" model for gp gp amplitudes 
were used, taking into account the nonperturbative proton wave functions. Even if we 
multiply the result of Ref. [12] by the suppressing factor, it will be larger than ours. This 
could serve as the indication of the role of nonperturbative effects. Our model is based 
on the Regge-eikonal approach for the amplitudes, which is primordially nonperturbative, 
normahzed to the data from HERA on 7p — > J/^p [6] and improved by the CDF data 
on the exclusive di-jet production [7] . 

To estimate the signal to QCD background ratio for bb signal we use the standard 
expression for gg — > bb amplitude and assumptions [13]- [16]: 

• possibility to separate final bb quark jets from gluon jets. If we cannot do it, it will 
increase the background by two orders of magnitude under the 50% efficiency. 

• suppression due to the absence of colour-octet bb final states 

• suppression of light fermion pair production, when Jz,tot = (see also [17], [18]) 

• cut Et > 50 GeV {9 > 60°), since the cross-section of EDD bb jet production 
strongly decreases with Et (see formulae (39)). 



The theoretical result of our numerical estimations is 

Signal{pp — > pHp — > pbbp) 



QCD background 



> 3.8 



GeV 
AM ' 



(24) 
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where AM is the mass resolution of the detector and Mh — 115 GeV, which can reach 
O.OlMjj due to application of the "missing mass method". Similar result was strictly 
obtained in [13], [14]. Under the above circumstances the total efficiency at the integrated 
luminocity 30 fb"^ is ~ 10% and numerically estimated significance of the event is about 
3(7, which is close to the one in 77 decay mode. 

5 Heavy quarkonium production 

Results for Xco.feo EDD production were obtained recently by some authors [13], [15], [19] 
in different approaches. To obtain the total cross-sections in the model considered in the 
present paper let us substitute widths of these states into (20). 

r(X6o ^ gg) - ^oixm ^ gg) (l + 9.8^) = 550 keV (see [13], [15] for details), (25) 

where the width is set to the lattice result ro(x6o ~^ gg) — 354 keV [20]. After replacement 
we obtain for LHC and TeVatron: 

(7pp_.p+;,,„+p ~ 1.3 n6 , = 14 TeV , ((3),(4) cuts), (26) 

(7pp_p+;,,„+p ~ 160 p6 , = 1.8 TeV , (CDF cuts), (27) 

The same procedure can be done for Xco- Taking the total width r(xco — gg) — 
14.9 MeV [21] we obtain 

~ 4 , = 14 Tey , ((3),(4) cuts), (28) 

(Tpp^p+^^o+p ~ 600 n6 , x/s = 1.8 TeV , (CDF cuts). (29) 

6 Radion production. 

Now there is a great interest to the multidimensional properties of the space-time. One 
of the models was proposed by Randall and Sundrum [22]. We have considered the case 
of one compact extra-dimension. In this case we have additional scalar particle Radion, 
that reflects the existence of an extra dimension and represents the fleld of the " distance" 
oscillations between the branes along the extra dimension. Since Radion has the same 
quantum numbers as the Higgs boson, they can mix [23]. After mixing we have two 
mass eigenstates, which could be observed experimentally. For the EDDE the following 
replacements in (22) should be done: 

fin) ^ 034/(77^*)+ 77 & for/i*, (30) 
/(??) ^ 7(«i2/(77.0 + 7a) forr*, (31) 

where 77^*,^* — ^ft* r*/^''^t other parameters are obtained from formulae in the Ap- 
pendix A of [23] : 
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7 = -y/A,^ , t; = 246 GeV is the Higgs VEV , A<^ is the radion VEV , (32) 
Z^^l- 6^7^(1 + 60 , tan 2^ = 12^7^-^— — 1 — ^ , (33) 

a^cose/Z ; 6 = - sin e/Z ; c = sin 9 - 6^7/Z ; d = cos ^ + 6^7/Z sin ^ , (34) 
ai2 = a + c/7 ; 034 = d + 67 , (35) 

2 22,22 2 j2 2|i,2 2 /'ocN 

rrij.* — c rrif^ + a rUj. ; nif^* — a rrif^ + rrij. , (36) 
r^ar* + bh* ; h^cr* + dh* . (37) 

Results for the total cross-sections in the case of Cgp = 3.3 with Sudakov-like suppres- 
sion are depicted in Figs. 5, 6 for several values of mixing parameter ^ and for the vacuum 
expectation value of the radion field — 1 TeV. 

7 EDD dijet production 

For the EDD production of a dijet system of mass Mx in the leading order (see, for 
example [4]) we have: 

I^PP-xr ^ E^4^ .^ = 99. (38) 

llirp 

\F,,-.x\' 3^^4 ,X = QQ,(3= ^1 - 4myMi , (39) 



and the cross-section becomes 



7 I 



Unit. 1 2 

(40) 



where = 4£;|/M|., and T^%jjp is calculated as in the section 2 with substitu- 
tions (38), (39). 

Note that all the results of this article are given for the value of the constant Cgp — 3.3, 
which is obtained from the upper bounds for the exclusive di-jet production at TeVatron 
energies [7] . Cross-sections and numerical estimations for Cgp at different transverse energy 
cuts are the following: 



Et> 7 GeV , a<3.7nb, Cgp < 3.3 (41) 
Et > 10 GeV , cr < 0.97 ± 0.065 (stat.) ± 0.272(sys.) nb , Cgp < 3.4 
Et > 25 GeV , cr < 34 ± 5 (stat.) ± 10{sys.) pb , Cgp < 4.2 . 

The lowest value is close to the result, obtained by fitting the HERA data on elastic J/^ 
production. It can serve as the indication of model applicability. 

From the analogous calculations for LHC with cuts (3), (4) we have: 
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Et > 10 GeV, 
Et > 25 GeV, 
Et > 50 GeV, 



a{pp — > p + jet + jet + p) ~ 7 n6 
a(j)p — > p + jet + jet + p) ~ 150 pb 
a(j)p — > p + jet + jet + p) ~ 8 pb. 



(42) 



8 Conclusions 

We see from the results that there is a real possibility to use advantages of the the EDDE 
for investigations at LHC. Accuracy of the mass measurements could be improved by 
applying the missing mass method [24]. 

The low value of the exclusive Higgs boson production cross-section obtained in this 
paper is mainly due to the Sudakov suppression factor (16), the full validity of which is not 
obvious, because the confinement effects can strongly modify the "real gluon emission". 
It is interesting that in spite of different models and quite different ways of account of 
absorbtive effects in our paper and in Ref. [8], the final results appeared to be quite close. 

Certainly, the cross-sections may be several times larger due to still not very well 
known non-perturbative factors. 

In the case of heavy quarkonium and dijet EDD production cross-sections are much 
larger than for Higgs boson production, and some other important investigations like mea- 
surements of the azimuthal angle dependence and the diffractive pattern of the interaction 
could be done. 
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Figure captions 

Fig. 1: The process p + p— + Absorbtion in the initial and final pp-channels 

is not shown. 

Fig. 2: The full unitarization of the process p + p— + 

Fig. 3: t-distribution da/dt/atot of the process p + p —>p + X + p for masses of the 
system X equal to 100 and 500 GeV. 

Fig. 4: ^-distribution da/d^/atot of the process p + p— ^p + X+ p for Mx = 100 GeV. 

Fig. 5: The total cross-section (in fb) of the process p + p^p + h* +p versus Higgs 
boson(/i*) mass for Cgp — 3.3 with Sudakov-like suppression at LHC. Parameters of 
RSI model are shown. 

Fig. 6: The total cross-section (in fb) of the process p + p —> p + r*+p versus Radion(r*) 
mass for Cgp — 3.3 with Sudakov-like suppression at LHC. Parameters of RSI model 
are shown. 
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Figure 6: 
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